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Abstract 

We consider the conductivity of excitations in short chains of optical cavities with two-level atoms 
in the models of JCH type, where either we explicitly take into account the photon jumps between 
atoms, or is merely a transfer of excitation from atom to atom. We found a non-trivial dependence 
of the conductivity on the intensity of runoff and inflow (quantum bottleneck) in the presence of 
dephasing noise (the effect of dephasing assisted transport). 
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1 Introduction and background 


Complex systems is the strongest challenge to quantnm theory since its inception. It is related to the 
fundamental problem of the complexity of the description of many-particle ensembles at the quantum 
level, which led to the idea of a quantum computer. The most critical element in our understanding of 
the physics of quantum computing is decoherence, which can be described within the concept of an open 
quantum system (see [T]). 

Kossakowsky-Lindblad master equation (see |2] and |3]), gives the simplest description of an open 
quantum system, where the influence of environment is expressed by permanent measurements of the 
system and is represented as Markov random process. Despite the fact that in many cases this influence 
can hardly be considered as Markov process, and there are non Markovian approaches (see.. Eg, |1]; the 
analytical analysis of interaction with invironment can be found in 0 ), the dynamics in terms of the 
master equation has the simplest form, it makes it easy to compare result with the unitary dynamics of 
the expansion of the system, where we include some part of the environment to the considered system. 

Coherence means the presence of off-diagonal elements of the density matrix; quantum effects are 
always associated with coherence, and the effects of information nature as quantum algorithms, moreover, 
require even the existence of entanglement. Influence of environment usually leads to the suppression of 
coherence and treated as noise in the practical use of quantum effects (see m)- However, in some cases, 
on the contrary, the influence of the environment can help to maintain the desired quantum effects. The 
most famous example - quantum Zeno effect, in which the measurements allow to conserve the desired 
quantum state (see the work jH] Zeno effect in NMR- quantum computing as the freezing of quantum 
states, coherence protection by Zeno effect in E)- It is also possible to suppress decoherence by a special 
noise (see, for example [8j). In the work [9]the non-trivial relationship between the unitary dynamics 
with Lindblad master equation is described, in particular, the possibility of noise smoothing by a unitary 
control. 

The influence of noise is very important for the conductivity of excitations in groups of two-level atoms 
enclosed in the optical cavities with the possibility of exchange of photons between the cavities. Excitation 
conductivity can be used to analyze the transmission of energy and information at the quantum level, in 
artificial devices and living organisms (see, for example [10] , jUj , [H]) where the addition of the noise 
to the model makes it more realistic. 
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We can formulate the conductivity of excitations in terms of quantum walks m, m), which have 
a valuable difference from classical random walks. For example, quantum conductivity quadratically 
exceeds classical: the conductivity of the classical walk on a linear chain takes time of the order of 
squared length of the chain, whereas quantum - of the order of the length ( HSI, [H], US]). For the walks 
on glued trees quantum speedup of the conductivity will be even exponential(see |16|). 

We will consider two important quantum effect associated with conductivity: dephasing assisted 
transport (DAT) and quantum bottle-neck. DAT plays an important role in the light harvesting Fenna- 
Matthews-Olson (FMO) complex in green sulfur bacteria (see |T7]). The essence of the DAT effect is 
that the presence of noise at the resonant frequency increases transfer of excitations, so that there is an 
optimal intensity of the noise when the conductivity is maximal. In the works |18|.|12|. |n],|l9],|20],|2l]) 
DAT effect was investigated on the model of excitation transfer, without the explicit consideration of the 
interaction between photons and atoms. 

However, artificial reproduction of excitation transfer is important, for example, for the effective 
conversion of light into chemical energy by optical elements. To simulate such systems we must take into 
account the interaction of atoms with light explicitly, as is done in this paper. 

Quantum bottle-neck effect is that the maximum conductivity occurs at the optimum value the force 
of interaction of the end node with the sink (flow rate). The excess of this value does not lead to an 
increase in conductivity (as in the classical case where the conductivity increases monotonically with the 
intensity of runoff), but decreases it, until the complete disappearance. This is like DAT, purely quantum, 
and, moreover, counter-intuitive effect. It is noted in the paper [20] (p. 10) that the bottle-neck can exist 
only in the case of coherent exchange between nodes, i.e. it has quantum nature, but there were no exact 
boundaries for Lindblad dephasing, in which this effect arises. In this paper we present such boundaries. 

The physics of quantum effects in the conductivity of excitations is extremely important in biology, 
as well as in quantum technology, so it is of interest to study these effects in more detail in the framework 
of Jaynes-Cummings-Hubbard model (JCH), introduced in , [^. The explicit consideration of 
photons and their interaction with atomic excitations explains many quantum effects important for the 
applications. These effects can be used in quantum computing: optical CNOT gates ([21]) , quantum 
computations on optical gates (see |25| . and |26| 1 - quantum computations on linear optical gates and 
measurements, photonic quantum memory (see, for example, |27|h the description of the other quantum 
effects and their possible applications can be found in [^ , [^ , [30] , |3T] , [32] , |33| , [31], [35] . 

Excitation conductivity has more detailed description in JCH model, which depends on its parameters: 
amplitude of photon jumps between cavities and amplitude of atom-photon interaction inside of one cavity. 
For example, in the work |36| excitation conductivity is considered in JCH model without phonons. We 
make such analysis for the model with phonons and establish the influence of atom-photon interaction to 
DAT effect. 

Our goal - the quantitative description of two quantum effects arising in transfer of energy along 
the chain of atoms concluded in optical cavities: quantum bottleneck and dephasing assisted transport 
(DAT). We formulate a model in which it is possible to combine the simultaneous influence of these effects 
on the excitation conductivity, and then present the results of numerical simulations. 

We also give a brief qualitative explanation of these two effects, and indicate possible ways to use 
these effects to improve the quantum models of optical systems associated with the conversion of light 
energy into chemical. 


2 Quantum effects in the energy transfer 

We give a definition, starting with the most detailed model "excitations, photons, phonons." We 
consider a linear chain consisting of identical optical cavities with two level atoms inside. Each cavity 
holds the photons with frequency lUc- Fock state with n photons inside i-th cavity we denote by |n) fi.. 
Inside each f-th cavity there is one two-level atom, which eigenstates: ground and excited we denote by 
|0)at I, |l)at i correspondingly. The difference between frequencies of these states (detuning) d = ojc — oJa 
is small in comparizon with each of them: d uj^ that allows to use rotating wave approximation for the 
atom-photon interaction at the large enough time frame, and Jaynes -Cummings model (JC). In addition. 
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each atom is placed in a bath of thermal phonons that have dephasing effect on its excitation. In this 
paper all phonons have the same frequency ujp, for which DAT effect is maximal. In fact, this frequency 
is equal to the difference between eigen frequencies of the Hamiltonian of JCH model (see below) that 
are responsible for the transfer of excitations (see the work | 12 | 1 . 

The transfer of energy from the cavity to cavity occurs through the flight of photons between the 
cavities where is the amplitude of the flight from the cavity j to the cavity i per unit of time, so 

that the amplitudes of the opposing flights are mutually conjugate, that is Sij = Sj^i. 

The Hamiltonian of our model is thus obtained by adding to the Jaynes-Cummings Hubbard Hamil¬ 
tonian Hjch the term of exciton-phonon interaction Hint ep' 

^ ^JCH T ^int epi 

Hjch = huta Y. Y + Yh^^ti 

i i i 

Hint ep = gYi^t + hi)afai. 
i 

Here af, at denote operators of creation and annihilation of a photon in Tth cavity, bf ,bf- denote creation 
and annihilation of a phonon in i-th cavity, denote creation and annihilation of the excited state 

of the corresponding atom (excitation), 7 is the amplitude of photon emitting by excited atom in the 
unit of time. We assume that the photon jump is possible between the neighboring cavities only where it 
occurs with the same amplitude so that all dij = 0 for |i — jj yf 1, and for |i — jj = 1, 6ij = S. Constant 
g is the intensity of exciton-phonon interactions (square root of Huang-Phys factors). 

In the simplified model, which we call ’’excitation-phonon” photons are ignored, and excitations are 
transmitted from one atom to the other; its Hamiltonian Hgp has the form 

Hep — He Hint ep-, He = hiOn ^ ^ ^i T A (^) 

i iiij 

The inflow of the energy can be viewed either as a constant occurrence of the excitation on the first 
node, taking the energy income from the outer bath, or simply set the initial state of the first atom as 
excited, regardless of photons. In the first case the bath inflows can create excitation irreversibly, or we 
can simply consider the initial state of a strongly excited field in the first cavity, which interacts with the 
first atom. 

There are two ways for the simulation of dephasing as well. Either we consider it as the interaction with 
the explicit phonons, or as the irreversible process with only excitations and without explicit phonons. 

Combining all of the above methods of studying the conductivity, we get all kinds of particular 
computing models. In any case, for the irreversible models master equation of Kossakowsky-Lindblad 
should be applied: 


i ■ St 

Pt+st = Ugi-ptUst + St''^{LipL* — -{L*Lip + pL*Li)), Ust=e ^ . (3) 

i 

which contains the unitary dynamics Ugt as the particular case. Here Lindblad operators Li describe 
the irreversible part of the process. For example, the irreversible runoff of excitations to the sink from the 
last node (end) is expressed by the operator Lgink = gsmfe| 0 )end|l)smfc( 0 |smfc(l|end, where the positive 
coefficient Psink expresses the intensity of runoff, dephasing in i-th node resulted from the interaction 
with the implicit phonons is reflected by the operator Ljeph, i = photon leak through the walls of 

the i-th cavity - by Ljet, i = Q-i, etc. 

We consider the resulting quantum effects of excitation conductivity. We determine the conductivity 
by the degree of filling sink in a certain time. 

1). Quantum bottleneck. Reduction of conductivity with increasing intensity of the runoff. Graph of 
conductivity depending on the rate of the runoff is not monotonically increasing as would be in the case of 
classical conductivity, but has a local maximum at a finite rate . The nature of this counter-intuitive 
effect is purely quantum. 
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Figure 1: Dynamics of the population of the initial state with a large runoff intensity. Runoff shifts the 
density matrix to the initial state, the population of excitation in the first atom (for two-atom chain this 
is Pin,in) is shifted to the maximum point, and remains almost constant for a long time, preventing the 
sink filling. 


It can be explained by considering the dynamics of the population of the initial state according to the 
time when the runoff intensity is large. If there is no runoff, the dynamics will be similar to the graph of 
the cosine, which has a maximum point at the initial time. The runoff constantly destroys the excitation 
at the end of the chain that makes the density matrix close to the initial state in which a reduction of 
the excitonic population at the first atom decreases with time very weak. This process competes with 
the obvious classic excitonic population decrease with increasing runoff. 

If the runoff is strong enough, it constantly holds the density matrix (and hence population of ex¬ 
citations in the first atom) in the initial state in which there is practically no reduction of excitonic 
population at all that results in ” freezing ” of the population, like in quantum Zeno effect, in which the 
effect of” freezing ” comes from frequent measurements of the quantum state. In the case of the classical 
conductivity instead of cosine fall of excitonic population we would have a linear decrease, and there 
would be no bottleneck. 

2). Dephasing assisted transport - DAT. A plot of the conductivity on the intensity of the noise (the 
number of phonons of the resonant frequency, or intensity of the interaction between excitations and 
phonons, expressed in coefficients of Lindblad for dephasing) has a local maximum at a non-zero point. 
This means that there is some non-zero noise intensity at which the conductivity is maximal. There 
are quantum processes, connected with the conductivity, the effectiveness of which is enhanced with the 
increase of noise, for example, mixing at random walks on graphs m)- Effect of DAT stands out among 
them due a special role it plays in biology. 

DAT as a purely quantum effect. Dephasing is direct suppression of the off-diagonal elements of the 
density matrix. Quantum dynamics is obtained by adding the amplitude states, the evolution of which 
leads to the same point of the classical space. Dephasing shifts the phase of the amplitude to a certain 
area. If this area would be narrower than the natural spread of the phase of states in the absence of 
dephasing, this influence will make the interference more constructive, thus increasing conductivity. This 
effect can be illustrated by the example of road traffic. If possibilities of cars are the same and allow 
rapid acceleration and braking, then synchronous mode of motion (phase conservation) will support high 
bandwidth of the road. But if the possibilities of vehicles are very different, then to increase bandwidth 
it is needed to the contrary, smooth mode of braking and acceleration, i.e. dephasing. 
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Qubit-based model 


3.1 Evolution equations 

Kossakowski-Lindlbad equation in diagonal form is used for calculating non-unitary evalution: 

= [H,p]+i -/,{AipA* - -{A*Aip + pA*Ai)) (4) 

i 


Unitary evalution is optimized by diagonalizing the Hamiltonian. 

In both cases, the evolution operator is determined by a Hamiltonian H with a set of non-normalized 
matrices Li = y/%Ai. 

For general (non-unitary) case, the evolution is computed by the equations ^ directy. 

These evolution equations, as well as the measurement implementation, are independent of the used 
model. Various models could be created by constructing a Hamiltonian H and a set of Lindblad matrices 
Li according to that model rules and semantics. 


3.2 Qubit-based model construction 

In order to allow simple operator-based constructing of Hamiltonians and Lindblad matrices we build 
the models by binding one a several qubits (a qubit group) to a specific state parameter, for example to the 
excitation of the first atom (see [3H]). The model construction methods work with virtually huge matrices: 
for example for a chain of one-levelled excitons of a length 100, it would have a virtual Hamiltonian of size 
(one qubit per each exciton). If excitations were four-leveled, that would make two qubits for 
each excitation. These matrices are not actually constructed in memory, an energy-limiting projection is 
used which selects only those states that have the energy between the specified lower and upper bounds. 
The resulting Hamiltonian and Lindblad matrices which are used in the actual calculation have much 
smaller sizes due to this projection. 


3.3 A chain of excitation-photon-phonon cavities 

A general chain model was constructed. This model implements a chain of cavities, each of which has 
simular multi-level excitation, photon field, and phonons. It has energy limiting support and supports 
attaching an input source and an output sink. 






Figure 2: The exciton-photon-phonon cavities chain schematics. 


There are three types of qubit groups in the chain: 
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• bi is the phonon state for atom i, 1 — the corresponding phonon is present, 0 — the corresponding 
phonon is absent; 

• Gi is the exciton state for atom i, 1 — excited state, 0 — not excited; 

• Pi is the photon state for atom i, 1 — there is a photon, 0 — no photon. 


The chain model has the following parameters: 

• Natoms — the number of atoms, 

• k — photon tunnelling rate between neighbouring cavities, 

• /i — photon-atom interaction strength, 

• g — phonon-atom interaction strength, 

• iOa — frequency of one excitation (atomic transition frequency), 

• LOp — frequency of one photon (cavity frequency), 

• Wg — frequency of one phonon, 

• in — input rate, replenishment coefficient for the first excitation ai, 

• out — output rate, sink coefficient for transferring excitation from to s. 

Hamiltonian of the system equals to Jaynes-Cummings-Hubbard Hamiltonian: 

H =Y^ u}pP+p~ + ujaaf a” + Y ^bbf K + Yi^Pi+iP^ + ^*Pt p7+i) + 

i i i i 

Y^pPi + p*pt0-7) + (g + 9*) Y^^^i + )ot0-7)■ 


(5) 


The sink s is attached to the last exciton in the chain. Input and output is performed using Lindblad 
operators: 


Lin = in *pf 
Lout = out * 

In all the follow-up graphs and results, it is taken Wq = ujp = 0.1 and ujg = 0.01. 


( 6 ) 


3.4 Dephasing models 

Two dephasing models were implemented: 

• Unitary-based: with explicit phonons as shown on|^ 

Hamiltonian of the system has the {g + g*) + ^t)ot0-7) part. 

i 

• Lindblad-based: there are no explicit photons in the system. Hamiltonian of the system does not 
have the corresponding part, Lindblad-like dephasing operators Di = g*pfp7 are used instead [TT]. 

With the addition of Lindblad-like dephasing operators, evaluation equation takes the following 
form: 

pt+st = U^tPtUst + StYiW - l{L*L,p + pL*L7)+StYiD*D,pD*D, - \{D*D,p + pD*D7) (7) 
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4 Quantum bottleneck 


Graph shows the optimal value of the input rate depending on the output rate for an chain of two 
cavities with one excitation and one possible photon each, for sample k and /r values. The second part of 
the graph shows the time taken to reach target sink value using the optimal input rate for the selected 
output rate. It could be seen that even with optimizing the input rate, the quantum bottleneck effect is 
observed. Graphs are capped at timeReach = 400, this is the reason why the first part of the graph is 
incomplete for ^ = 0.1. 

Graph shows the optimal value of the output rate depending on the input rate for the same chain 
of two cavities. 


inOpt 



The optimal input rate, over the output rate 


timeReach 



Time taken to reach target sink value 0.995 
using the optimal input rate, over the output rate 


Figure 3: Optimal input rate for for a chain of two atoms. 


outOpt 



The optimal output rate, over the input rate 


timeReach 



^=0.1 k=0.1 
^=QA k=0.5 
IJ=Q.B k=0.5 
^;=0.5 k=0.8 
^=0.8 k=0.8 
^J = ^.Q k = 1.0 


Time taken to reach target sink value 0.995 
using the optimal output rate, over the input rate 


Figure 4: Optimal output rate for for a chain of two atoms. 



Figure 5: The dependency of time taken to reach target sink value 0.995 over in and out rates, /r = 0.8, 
k = 0.5 
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outOpt 


0.2 0.4 0.6 0.8 1.0 

The optimal output rate, over ^ 




- ^J = OA 

- f/ = 0.2 

- f/ = 0.5 

- IJ = O.S 


The optimal output rate, over k 


— k = 0.1 

— k = 0.2 

— k = 0.5 

— k = 0.8 


Time taken to reach target sink value 0.995 
using the optimal output rate, over /r 


— fj =0A 

— IJ =0.2 

— IJ =0.5 

— fj = 0.8 



Time taken to reach target sink value 0.995 
using the optimal output rate, over k 


Figure 6: Optimal output rate for for a chain of two atoms without input (initial state has a photon in 
the first cavity). 






Figure 7: Evolution of the sink state over time 


Note that all two-dimensional timeReach graphs in the paper have logarithmic scale. 

For all numeric experiments with input enabled the total energy was not explicitly limited, except for 
the natural limits imposed by the model. For numeric experiments with input disabled, the initial state 
has one photon and the total energy is limited to 1 (because it could not raise higher than that). 


































On figure]^ it is seen that for a chain of 2 atoms the optimal output rate has little dependency on the 
photon tunneling rate, but high dependency on photon-atom interaction strength. 

It is also visible that increasing the photon-atom interaction strength over some critical value could 
have a negative effect on conductivity. 

From the evolution graphs in figure it could be seen that the optimal input/output rate depends 
on the target sink value (though that dependency is not critical of high enough target sink values). For 
example, with low target sink value 0.3 and k = 1.0, /i = 1.0, it is more effective to use in = 1.5, out = 2.5 
than in = 1.5, out = 1.5. 

5 Dephasing-assisted transport 

5.1 Lindbad-based dephasing model 
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fi = 0.2, k = 0.8, time = 150 /i = 0.8, k = 0.2, time = 60 

Figure 8: The state of sink at the fixed time. Two atoms, no input (initial state has a photon in the 
first cavity). Dependency over the output rate out and the dephasing coefficient g. 

For high enough times, the absolute maximum over all possible output rates is reached when g = 0. 
But non-zero values of g can greatly improve the conductivity in cases when the output rate does not 
match with the optimal output rate. That is shown on figure It could also be seen that the effect does 
not depend on g being equal to some exact value in some cases, and there is a broad range of possible 
values for g that improves the conductivity. 

For very low times (or low target sink values, depending on the stop criterion) these graphs are 
different and could show other short-lived effects, but we are inspecting target sink values close to 1.0. 

In the first part of (/j, = 0.2, k = 0.8) it could be seen that non-zero values for g improve the 
conductivity for the case out < outOpt. In the second part (/i = 0.8, k = 0.2) it could be seen that 
non-zero values for g improve the conductivity for the case out > outOpt. 

The exact range of possible output rate values for which the conductivity could be improved by adding 
dephasing depends on the parameters of the chain. 

The same effects could be observed for higher chain lengths, as could be seen on figure El 
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sink 


sink 



Figure 9: Evolution of the sink state over time. Two atoms, no input (initial state has a photon in the 
first cavity), k = 0.2, fi = 0.8. 



Figure 10: Evolution of the sink state over time. Two atoms, no input (initial state has a photon in the 
first cavity), k = 0.8, ^ = 0.2. 
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Figure 11: Evolution of the sink state over time. Five atoms, no input (initial state has a photon in the 
first cavity), k = 0.8, fi = 0.2. 


5.2 Unitary dephasing model 



Figure 12: The state of sink at the fixed time. Two atoms, no input (initial state has a photon in the 
first cavity). Dependency over the output rate out and the dephasing coefficient g. 


For the unitary dephasing model, the main observed effects are the same as for the Lindblad-based 
dephasing model in |5.1| — non-zero values of g can improve the conductivity in cases of inoptimal values of 
the output rate, though the overall effect of unitary dephasing (with the supplied parameters) is milder. 

Figure [T^ depicts the same experiment as figure but with the unitary dephasing model. 

On figures [^and it could be seen that setups that were optimal to reach low target sink values 
(for example 0.3 or 0.4) are not always optimal to reach target sink values close to 1. 
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Figure 13: Evolution of the sink state over time. Two atoms, no input (initial state has a photon in the 
first cavity), k = 0.2, ^ = 0.8. 






Figure 14: Evolution of the sink state over time. Two atoms, no input (initial state has a photon in the 
first cavity), k = 0.8, n = 0.2. 
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6 Conclusions 


We have shown the efficiency of the qubit representation of a JCH model with the Lindblad-like form 
of dephasing, and its relation with the simple unitary model of dephasing. 

Quantum bottleneck effect was explained and reproduced for both input and output rates. It was 
shown that, depending on the model layout, the optimal input rate could mismatch the optimal output 
rate. The dependency of the bottleneck effect on the in-model interaction rates was shown. 

Two models of dephasing showed similar results in the long term, when the numeric experiments were 
run long enough for the sink to reach a value close to 1. For both models non-zero dephasing did not give 
a positive effect on the conductivity in the case of optimal output rate (when the quantum bottleneck 
effect is not visible, but gave a large positive effect in some cases when the output rate was not optimal 
and when the conductivity is capped by the quantum bottleneck effect. 

It was also shown that short-term effects differ from the long-term effects, and while some settings 
are optimal for reaching a low target sink value (for example 0.3), they could be far from optimal for 
reaching a close to 1 target value. 
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